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CHAPTER I 


SUMMARY 

A multifluids model to investigate ionospheric dynamics has 
been established on kinetic theory. Its resultant equations are used 
to examine the following dynamic problems in the y region of 80- 2000 Km 
of the ionosphere; 

i) Propagation of acoustic modes in the 500-2,000 Km of the 
ionosphere (two fluid model), 

ii) The relation between the cross field plasma drift 
instabilities and type I and type II ionospheric 
irregularities. 

iii) Time dependent neutral wind structure and horizontal 
pressure gradient. 

In these studies, we have found out that case (i), the dissipation 
of acoustic waves by electrons, is through the viscous and thermal con- 
duction effects and ion is through Landau damping in the region of 500- 
2,000 Km of the ionosphere; the type I ionospheric irregularity (case(ii)) 
can be excited by either temperature gradient or density gradient along 
the vertical height, and the type II ionospheric irregularity can only 
be generated by temperature, gradients; case (iii), the ion drift, 
will initiate neutral motjof!- Ir an early stage of development; how- 
ever, the neutrals will also drive the ions in a later stage. 



CHAPTER II 

TWO-FLUIDS DESCRIPTION OF PROPAGATION OF ACOUSTIC MODES 
IN THE REGION OF 500-2000 KM IONOSPHERE* 

II-l General Statement of the Problem 

Propagation of acoustic waves in the upper atmosphere of the earth reflects the 
thermal fluctuations of the motions of the upper atmosphere. The investigation of the 
characteristics of the propagation and dissipation of acoustic-mode- waves Is of 
particular interest to environmental perturbation and to plasma physics laboratory 
definition and also in the design and operation of the space shuttle. 

The physical properties of the medium in which waves propagate affect the 
characteristics of the wave propagation. In particular, wave propagation is 
governed by the transport phenomena of the medium which depend on collision 
effects of particles and waves. As the acoustic modes propagate in the various 
altitudes of the atmosphere, the characteristics of the wave propagation vary from 
region to region because the physical parameters of the medium also change. In 
the lower altitude region, the atmosphere contains more neutral molecules than 
ionized particles. As altitude increases, the particle ionization increases. In the 
meonwhile, neutral -neutral, ion-neutral and electron-neutral collision frequencies 
decrease. As the altitude reaches 500 Km, the electron collision frequency 
is400sec \ the ioncollislon frequency is 0. 1 sec"*, and theneutral collision frequency 
—2 —1 

4 X 10 sec (Hanson 1965; U.S, Standard Atmosphere, 1962), In other words, 
the effect of neutral particles on the propagation of acoustic waves with wave 
periods 0,01 to 10 seconds (or wave frequencies 0,1 to 100 Hz) is negligible 
in the altitude above 500 Km. Thus, at altitudes above 500 Km, the collision 
effect due to neutral particles disappears, and we have to consider Coulomb 
collisions instead of neutral particle collisions. 

As we have mentioned earlier, propagation of waves is modified by transport 
effects due to either particle collisions or Coulomb collisions. In the lower’altitude 
atmosphere where the neutral particles dominate, the transport coefficients of 


*A portion of the results are published in the J. of Geophysical Research 
(1975). 
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viscosity and thermal conduction from kinetic theory (Chapman ond Cowling, I960) 
are 


K c^m(9y-5) 

where n, K, n, and T denote viscosity, thermal conductivity, number density 
and temperature, respectively, and C^, y and denote specific heat at 

constant volume, the ratio of specific heats, and the neutral-neutral collision fre- 
quency, respectively. In the collision-dominated regions, such as the D, E, and 
lower F regions of the ionosphere, order of magnitude estimations show that both 
thermal conduction and viscosity have a negligible effect on the profxjgation of 
waves with periods from O.OI to 10 seconds. As the altitude increases, the 
collision frequency becomes smaller and both thermal conduction and viscosity 
become significant. Thus, at the higher altitudes the acoustic modes suffer from 
strong dissipation effects due to thermal conduction and viscosity. This Implies 
that the acoustic modes might not propagate very long distances without being 
damped out. It seems that this conclusion was reached too fast. In 
order to have a better understanding of the physics of acoustic mode propagation 
In the upper atmosphere, particularly in the transitional region from the collislohal 
to the collisionless situation, a new approach based on kinetic theory should be 
considered. 

As our primary interest is in investigating the acoustic modes propagation 
In the upper atmosphere where particle transition from collisional to collisionless 
conditions occurs, our problem is to Investigate modes with wave periods from 
0.01 to 10 seconds propagating In the altitude region from 500 Km to 2000 Km . 

In this region, isotropic electrons are in transitional between the collisional to 
coillsionless states, while the anistropic ions are In a collisionless state. Our 
primary consideration shall be limited to wave propagation in the medium of the 
ionized plasma rather than the neutral particles because the effect of neutral 
particle collisions is negligible in this region. The propagation of acoustic modes 
in the ionized particles has been studied for several years. For the case of wave 
propagation /n a collision-dominated fully ionized plasma, a comprehensive study 
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was done by Braglnskfl (1965). For fhe acoustic mode propagation In a collisionless 
fully ionized plasma, interesting results were obtained by Bernstein and Trehan 

(I960), In a study of acoustic wave propagation In the medium of ionized particles 
without a magnetic field, Bhatnagor, Gross and Krook (1954) presented a model of 

single species rarefied ionized particles, and Gross and Krook (1956) showed a 
model of two-component ionized gas. In the study of the dynamics of rarefied neutral 
particles, Hamel (1965) accomplished an investigation of a binary gas mixture and 
Holway (1966) proposed a statistical model with a multicomponent gas mixture. To 
investigate the acoustic wave propagation in the transitional ionosphere, Tolstoy 
(1972) suggested that continuum theory is inadequate and a more sophisticated 
treatment based on an adequate molecular model combined with Boltzmann's 
equation should be adopted. 

In the present study, we propose to investigate acoustic mode propagation 
in the upper atmosphere from 500 to 2000 Km altitude, where electrons behave as If 
they are in a transition state between a colUsional and a collisionless state, 
while the ions act as though they are in a collisionless regime for waves with 
periods of 0.01 to 10 seconds. The rigorous approach io the analysis of 
transitional electrons from the case of collisional to collisionless by using kinetic 
equation is almost impossible. However, it has been shown by Hung and Barnes 
(1973b) that a mathematical similarity exists between the electron viscous damp- 
ing and Landau damping as the electrons transition from a collision to a colli- 
sionless regime. *In other words, the electron viscous damping changes con- 
tinuously into Landau damping as the ratio of wave frequency to electron colli- 
sion frequency passes through unity. A similar conclusion was given by Rognlein 
and Self (1971) based on the hydrodynomic description. Therefore, in the present 
analysis, the dynamical behavior for the transitional elec;'rons can be described 
by the fluid-like equations, while the anisotropic ions follow kinetic equations 
modified by ion-electron collisions. 

Before closing this section, let us discuss the validity of the present theory 
in this transitionol region. Assuming that ion and electron temperatures are in the 
same order, the criteria for the validity of the present theory is 



where (1), T and t. denote wave frequency, electron and ion collision time, 


0 ) 


respocUvely. There arises a quesHon of whether the fluid-like equations are 
proper to describe the behavior of electrons under the criteria of validity 
mentioned in Equation 0 )* To answer this question# lot us look at the conditions 
for electrons to bo able to satisfy the fluid equation* The criteria for electrons 
to behave like a fluid shall satisfy the following conditions: 

« I i I (2a) 


jki. 1 a T « 1 
e e 



where a and a. denote electron and ion thermal speed# respectively# and the 

0 I 

subscript H implies the component along the geomagnetic field. In the wave- 
particle interaction problem# it Is generally true that utA o and uVkS. a. . 
By using this choracteristic and the condition (2a)# ^ 1 we may have 

Ikjila^ ~|u) land jj 1 <*^1 when waves propagate fairly near 

'^/ 2 * Thus# this deduction gives the conclusion which is essentially the 

criteria (2b). In the limit jiol t 1# as pointed out by Barnes (private 

communications# 1974)# the present theory at least gives qualitatively correct 

results when luJ 1 t a few tenths, and the coilisionless theory should be 
e 

qualitatively correct when | u) jt 5 3. Therefore# the similarity between 
viscous and Landau damping suggests a continuous transition between the two. 

In this sense# the colllsional term is proportional to out # and the Landau damping 

<-•1 ^ • 
term is porportional tom (k , a )** • Thus# the two terms are always comparable 


■|he ditterence in behavior between collisional and collisfonless electrons lies 
in (I) isotropy of electron pressure tensor# (2) the heat equation# (3) the 
frictional force# and (4) viscosity . Bornes(private communications# 1974) 
indicated that only the collision time# not the mean free path# is relevant for 
collisional isolropization# heat conduction along the field lines, which has a 
scale 1A|| (not lA), is the only appreciable heat transport process; the 

frictional force (in this case '*^ ^|| T^) is governed by the length scale lAjj • 

Thuslsotropy# the heat equation# and friction are covered by fluid equations 
If Equation (2b) Is satisfied. Viscosity Is more complicated# but the con- 
clusion is the same . 
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at the time scale of transitlonol from collislonol to collisioniess electrons. 

As the altitude decreoses below the present interest/ olectron-neutrol 
and ion-neutrol collisions become significant/ and wo have to treat as a three 
fluid model (e.g . Hung and Wu,l97/i.),v 
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II“2 Mathematical Formulations 

To Idealize the problem In hand, assume that the plasma and magnetic 
field are, on the overage, uniform throughout an effective!/ Infinite volume, that 
there is no net average electric current, and that the average electron pressure 
tensor, but not the ion pressure tensor, is isotropic. Lot jv, P, jg^and n denote 
velocit/, pressure, heat flux, and viscous stress tensor respectively and lot sub“ 
scripts or superscripts of e and I denote electrons and ions for the proper situation, 
respectively. £and JB represent the electric and magnetic fields, e the Ion charge, 
m, the Ion or electron mass and c the speed of light. Then the fluid equations 
fo^ the electrons and the ion velocity distribution f. may be written as follows 
(Braginskii, 1965; Hung etal., 1973): 

(3) 

(4) 



m n 


dv 

_ 


e'e dt 


- 7 p - 7. n® - en (E + — 
^ e w e c 


v^ X 


R 


and 



dT 

e 

T 


+ p 7 
e ~ 



(5) 


3f. 

i_ 

St 


+ V.7 f. + — (E+— 'vx B) 

.s,. I m ^ r 


m. 

I 



m 

e 


m.T 
I e 








1 


m. n. 

I I 


Sf. 





where P *n T*, p/, = n. is the electron 

e ® ® (fi) * ~e ~ e 

collision time, and R “ R. is the coHisional momentuim transfer from ions to 
~e ~i 

* 

Here the unit of T is in erg 
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electron and vise versa . £ Is composed of frictional force and thermal force 
R T/ R ” £u electron thermal fluxjg^ Is composed of onalogous parts# 

= q^+ , Finally# of course# the electromagnetic fields must satisfy Maxwell's 

equations. 

As it stands, Equations (3) - (5) show the governing equations of electrons# 
and Equation (6) is the governing equation of tons. Equation (4) Includes the 
colllslonal momentum transfer from ions to electrons# and we have Ignored colllslonal 
energy transfer from ions to electrons in Equation (^. This Is because the colllslonal 
ion“electron energy exchange, whose characteristic time is on the order of or greater 

1/2 

thanCmj/me) t. which Is greater than the time scale of present interest# 
has been neglected. The rlght-^hfnd side of equation (6) shows the ion- 
electron collision term which Is of the same form as the Fokker-Planck colllslonal 
term that describes Brownian motion of particles In a moving medium with temperature 
T . T(>e first term of the collision terms describes the colllslonal energy transfer from 
electrons to Ions# and the second term Implies the colllslonal momentum transfer from 
electrons to ions. Again w« are going to neglect colllslonal energy transfer term 
because the time scale of the colllslonal energy transfer which affects the evolution 

I /2 

of distribution function is on the orderof orgreater than T (mj/me)or t. (mj 
that is long compared with the time scale of the problem of current interest. 

In the present analysis, we are concerned with waves whose circular fre- 
quency lU Is small compared with the Ion gyrofrequency# and whose wavelengths are 
long compared with the mean Ion Larmor radius. Under such circumstances# the 
momentum transfer due to colllsloni can be reduced to (Cf. Hung and Barnes, 1973a; 
1973b; 1973c; Hung et ol# 1973 ) 


R = -0.71 n K| T 
-«e enl e • 


( 7 ) 


where the subscripts II# X refer to the magnetic field direction e ~ B/ | B | . 
Similarly# the electron heat flux Is 




( 8 ) 
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whoro is coefficlonf of elocfron heat conductivity. Furthermore/ under the 
' ' 0 

present condition/ the stress tensor becomes (Hung etol/ 1973) 


= ~ 2 


■ 

VO L 


r av° 

J_| 

(av® av® av®\l 

z 

>L + L + ) 

_ a z 

3 * 

^ ax ay 3 z / . 


n ®=TT °=-T]° 

XX yy VO 


Lax 


tt ® r-n ® 

xy yx xz zx yz 


\ax 


e - 0 . e 

X ^ ^y ° ^ 7 . 


ay a z 
(10) 




v^here = e B / m c is the electron gyrofrequency/ and under the current 
e o 

condition 0 t >> The zeroth order coefficient of electron viscosity Is 
e e 


11 *"= 0.73 n T T 
VO e e e 


For the ion kinetic equation/ the first term on the right- hand-side of Eq. 
(6) affecting the evolutic«n of the velocity distribution on the time scale Is.^^ 

1/2 

T (m. /m )' ' which is long compared to time scales of Importance for 

© I 6 

the wave/ and may therefore be neglected. Hence/ from Eq^ r/)/ the right-hand 
side of Eq. (6) is just (Hung et al/ 1973) 


m. n. 

I I 


li. order to solve Eqs. (3) - (6)/ we assume T - <T>+ ® T ^/ t)/ n “^n^ 
6n fc,/ t) / P= <P> -f' 6P (x/ t)/ B = <B>+ 6 B (x / t)/ E = 6 E (x , t)/ 


*In the ionosphere/ the situation is more favorable to this condition as the 
altitude is higher. In the present case/ let us choase the lowest limit of 
altitude/ ^y 500 Km/ ~ lO^ Hz/ and ~ 2.5xl0”3. Then/ 

,Q Te 10^ which is much greater thon unity . ThuS/ i t can be sure 
that in the current region of interest/ the condition t is alwoys much 
greater than unity. ® 



~ 6 V** (x / t)and f. ~ < f. > + ^ f, ( v / 0^ whore < > denotes ensemble 

rv ' I I l’^ 

average. Wo consider the limit of o small amplitude pcrturbafiotv l6n/<n> j« 1/ 
etc. If the fluctuations are sinusoidal/ the technique of Fourier-Laplaco transfor- 
mation is applicable. Following the calculation procedure of Hung and Games (1973o/b) 

wo will bo able to obtain velocity fluctuations of ions and electrons in terms of the 
mobility tensor M, 

' M j 


6 


V . 


c 

<B> 



6E 


j “ 0 or I 


(13) 


For the limit of hydromagnetic waves, the displacement current is always; negligible. 
Then the electric current density may be written as 


I 


jjM 

4 tt 


K 


6 E 


(M) 


where the dielectric tensor is 


e - 

w 


e < n > (M. - M ) 

B > « I « e 


(15) 


Substituting the mobility tensor M. from Eq. (13) in Eq. (15), the explicit form 
of the dielectric tensor is 


K 


XX 






K 

K 

yy 

yz 

K 

K 

zy 

zz 


( 14 ) 


where the components of the dielectric tensor are 


1 


XX 


K = 

yy 



2y: 

k 

k 


2 

* 

1 

2 

1 

2 


- 1 


1 


] 


el 


2y? e 


{ • 


2 2 


k^a"»)po‘C / 2(3r+i2C) 

+ i 4- ° ( 1 + « ^ 

’ 2 2 V a 

2 (1)^ ' 

e 

.2 


) 


zz 


» * Z'. . 2 

-0-A) -4-^ ^ -2f2 


I 2 
'll ‘‘11 


o k,, a 
e II 


' k. 0J_,^ kj 3r^ + i2C^ 

w .. U)0. a 

*a ' ‘•''H i i 


/pj \ k, 0) Z'. ( 1 

K ■ - j (-^ ) ‘ — 

* n 2k„ <«n^. 


«> k, 

31 J- 


“"i ‘‘11 




, r = l + 2/(3+2iTi ) , 


**' y„) 2'<ii'""i 

Here 

n^C=l<n K^/ ('“ <%>)■ a= 3(1.71 r^- 0.71)-4i , 

C = 0.73 UJ T , a,\. = 2 <T;i . >/m. , a^ = 2<T >/m / 

V e (j|) (j|) I e e 

pj , = 8n <n> < T', 1 ,>/< U = u)/lk,. | C. , k = {k,, 0, k ) 

(fr) , II A - J. M 

= < B > / (4TTfn. <n>) Alfven velocity/ yj ” ' |k|j I ^ / 

A = 2*13(r -l)/a / and Z. = Z 21. are the plasma dispersion 

e III 

function and its first derivotive (Fried and Conte/ 1961) respectively . 
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The requirement that Eq. (14) bo consistent with Maxwell's equations gives 
the dispersion relation 



Substituting Eq. (16) in Eq. (17) gives 




K 


2 2 

1 ' 


k I c 


zz 


lU 


- K K 1 = 0 


(18) 


to lowest order in u)/n. . The first factor is the usual Alfven wave dist 7 ersion re- 


lation 


fe) 


(19) 


and the second factor, which gives the magneto-acoustic dispersion relation, is 


{ 


1,2 2 

K - - “ K . K. = 0 


yy 


tu 




zz zy yz 


(20) 


since Ik I » 1 k c/i» . 


zz 


By using Eq. (16), we can rearrange Eq. (20) by straightforward calculation 


to gi ve 


j^i+ -j- (B^“ P,'i ) cot^e = s(u,(UT^) 


(21) 


where 


i \2 


-S(U,^rJ= 

Pi I 


P P, 

e 


1 


ft ® 

^11 

4 \ -.1 J. n JT> 2 


P C 

z' - I 




, C )z.' +2 (r+i-^c 

3 V I e 3 V 


If 1 

^ ■ ^11 * iJ 


4 a 


Mr - i C ) z. - 2 

e e 3 V I tl 


( 22 ) 


2 

and 9 = 2f. ( k , < B > ) , B = 8n <nxT >/<B> . 
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1 1-3 Special Cases of Dispersion Relation 

EquaHon (19) gives Hio dispersion relation for Alfven waves in o plasma with 
anisotropic ions. The magneto-acoustic dispersion relation, t .e. Eq. (21), is tran- 
scendental and must generally be solved numerically. In order to look at the physics 
prescribed, it is convenient to consider a few interesting special limits that can pro- 
vide insight into the character of the v^oves without laborious calculations. 

A. Propagation Parallel to ^ ^ 

When Jk 1 1 < cot^ 6 -•« , go that the left-hand side of Eq. (21) be- 
comes zero, i .e.. 



(?i- ^i, ) 


or S -• » (U finite), we find from Eq. (22), 




2<T' > 

<T > 
e 


r - 1 

e 3 V 


(23) 


(24) 


Equation (23) shows that one of the magneto-acoustic modes propagates as an Alfven 
wave when 0 = 0 in a plasma with anisotropic ions. Eq . (24) corresponds to a 
compressive wave propagating along ^ ^ ^ • The analogous wave in magneto- 
hydrodynamics is the sound wave, and in a collisionless plasma the analogous wave 

is the ion acoustic wave which generally undergoes Landau damping (Fried and 

« 

Conte, 1961). 

In the transitional case, the electron collision frequency is on the same order 
of magnitude as the wave frequency, is the order of unity . From the result given 
by Hung and Barnes (1973b), we found 

|t)^^ 1>> 1 and «1 (25) 
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Under this condition, the dispersion relation con be reduced to (Hung et al 1973) 


(t) 


- C 

m. s 

I 


where C denotes Ion acoustic velocity, and 


Imuj fa -0.49 (lUT ) C k 

6 S 


or 


Im.UJ fa*0.67 C C k 

V 5 


(26) 


(27a) 

(27b) 


Similar to Hung and Barnes (1973b), this result indicates that all ion acoustic 
waves propagating parallel to the geomagnetic field are Landau-damped, and the 
damping is weak only if <Tg> » Particularly, for the case in ionosphere, 

< T > is on the order of <T.> which implies that most of the ion acoustic mode 
shall be damped out due to strong Landau damping. Furthermore, Equation (27)also shows 
that the dissipation rate is proportional to k under the condition when ujTg is on 
the order of unity. This result is similar to the collisionless case obtained by Fried 
and Gould (1961). 

For the collisional case, let us consider when uiTg is the order of (m / m.)* 

e t 

The dispersion relation for the weakly damped solution becomes 



25 + 4ti_^^ <T > 
TC e 

. 2 m. 

15 + 4n^c 


Imu) _ 
ReO) 




(23) 


(29) 


To be consistent with the collisional assumption, Tly^ shall be small compared with 
unity. Under this assumption. Equations (28) and (29) become 


<T > 
e 


/ it' \ ^ 5 e 

\k) ” 2 m. 


3 S 


2 2 

ImOJfa -0.21 T a k 
e e 


(30) 

(31) 
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which is a sonic wove propagating parallel to <B > • Equation (31) shows that the 

2 ~ 

dissipation rate is proportional to k under ihc condition when u) T is on the order 

of (m^/m.). In other words, the dissipationo! rate obtained from the collisionol 

assumption is equivalent to the collision-dominated case (for example, see Braginskli, 

1965). Strictly speaking, under the condition, whenWT^ is on the order of 

(m /m.), ions are in the transitional from collisionol to collisionlcss, and the 
e 1 

arguments similar to Hung and Barnes (1973c) should be employed. 

Equations (27) and (31) show that dissipation rate of acoustic modes gradually 

2 

change from the factor proportional to k to k as the dissipation mechanism transit 
from Coulomb collisions to Landau damping. 

B. Propagation Transverse to < ^ 

When kX’^B>, Z! “*0 and rL_ -♦ 0. Under this condition, F 
becomes 5/3 and expression (21) simplifies considerably. 



(32) 


Therefore, we have 




<T > 
e 


m. 

I 


C <T > 

V e_ 

3 m. 


(33) 


In the transitional case, is the order of unity. Equation (33) can be 
further reduced under the assumption of weakly damped solution, i .e., 





2<T> 




<T > 
e 


m. 

i 
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or 


I m oj ~ “0»0888 







k2 


(35b) 


For the collisional case, C is small as U) t is much less than unity. The 

V e ' 

dissipation rate becomes 


Imi«~ “ 0.122 T 
s e 


(36) 


It Is interesting to note that, under present calculation, the phase velocity 
of the transitional case in Eq. (34) is intermediate between the analogous expression 
which arises when ions are collisional 



<T.>+ <T > 

I ^ 

m. 


and that which arises when electrons are collisionless 


(- 


(U 

Fy 


C 2 + 2 — ^ ± — 

A m. 


(37) 


(3B) 


In the meanwhile, the dissipation rate also changes from the factor proportional 

2 

to k to k as the wave particle interaction changes from collisional to transitional 
cases. In a collfsionless plasma, strong Landau damping can occur for 0 slightly 
different from n/2 (Barnes, 1966). However, when j< is exactly perpendicular to 
<B>, this Landau damping vanishes. In the transitional case, the result is different. 
No matter how small the electron viscosity is, as long as is not zero, viscous 
dissipation always can be found even when k is exactly perpendicular to <B> . 
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II-4 Conclusions 


In the present study, we have Investigated the mechanism of acoustic 
mode propagation and the characteristics of wave dissipation in the region 
of the upper atmosphere in which the transport phenomena are transitional 
from Coulomb collisions to Landau damping. The results from the current 
analysis are helpful to understand how the acoustic modes propagate and 
dissipate in the upper atmosphere from ionosphere to magnetosphere, which 
is significant in the environmental perturbation studies. 

As shown in Tables I and II, the present study of acoustic-modes 
propagation, based on a theoretical model of the upper atmosphere from 
500 to 2000 Km altitude, shows that the dissipation rate changes from a 
factor proportional to 'ihe square of the wave vector to a factor only pro- 
portional to the wave vector as the altitude becomes higher. These results, 
approaching the collisional-limited case, are similar to the collision- 
dominated case {Brapinskii, 1965) in which the dissipation rate is proportional 
to the square of the wave vector, while the limit approaching the collisionless 
case is similar to the collisionless case (Fried and Gould, 1961) in which 
the dissipation rate is proportional to the wave vector. Thus, the 
dissipation rate changes gradually from a factor proportional to the square 
of the wave vector to the wave vector as the wave-particle collisions change 
from collision-dominated to collisionless situations. Physically, the 
dissipation mechanisms are also different. In the current transition model, 
acoustic wave dissipation by electrons is through viscous and thermal 
conduction effects, and by ion particles through Landau damping (wave- 
particle resonance interactions); while in the collision-dominated case, 
it is simply due to viscous and thermal conduction effects, and in the 


colHs1on-free case 1t is due to Landau damping. For the waves propagating 
along the magnetic field, two extreme examples are given for the transitional 
model, namely the colllslonal-llmited case with on the order of (mg/m^) 
and the transitional -limited case with ojt on the order of unity. The 
colllslonal-llmited case of the transitional model shows that the square of 
the characteristic phase velocity of acoustic waves In square of Ion acoustic 
velocity multiplied by 5/3 or y {where y Is the ratio of constant pru^sure 
specific heat to constant volume specific heat) which looks quite similar 
to the collision-dominated model in which square of the phase velocity Is 
the square of ordinary sound wave speed, yP/p (where p = m^Og + m^n^ + 

In the transitional-limited case of transitional model phase velocity Is 
equivalent to ion acoustic velocity which Is exactly corresponding to the 
phase velocity in collisionless case. In other words, change of the phase 
velocity of acoustic modes from collislonal to collisionless cases are from 

C 

a factor of y or ^ to a factor of unity. Thus, it is expected that the 
phase velocity of acoustic wave decreases gradually from collislonal to 
collisionless cases if the temperature keeps a constant value. Similar 
results are also shown for v/aves propagating transverse to magnetic field. 

Finally, we may conclude from these studies that the long wavelength 
acoustic mode 1s a much less damped wave for the collislonal than the 
coll Isionless cases, and the short wavelength acoustic mode is more easily 
observable in the collisionless than in the collislonal cases since the 
dissipation rate of the acoustic mode changes from a factor proportional to 
the square of the wave vector to the wave vector as the altitude Increases 
in the transitional region of the ionosphere. In the meantime, the physical 
mechanisms and phase velocities are also different for collision-dominated 
transitional and collisionless cases. 
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Table I. Comparisons of the Propagation of Acoustic Modes in CollisTon*Dominated, Transitional 

and Collisionless Media Along Magnetic Field 
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CHAPTER III 

THREE FLUIDS DESCRIPTION OF THE CROSS FIELD PLASMA 
DRIFT INSTABILITIES AND EQUATORIAL SPORADIC E* 

III-l General Statement of fche Problem 

Triggering mechanisms of sporadic E, in general, depend on the 
geomagnetic latitudinal location. For example, sporadic E in the high 
latitude region and the impinging of energetic particles through magneto- 
spheric cleft, is considered to be responsible for the mechanism. In 
the middle latitude, it is believed that the redistribution of ionized 
particles caused by the wind sheur may be responsible for the creation 
of sporadic E. Concerning the equatorial region, Matsushita (1951), 
suggested that the equatorial sporadic E (Es-q) can be related to the 
equatorial electrojet. Since then, intense study of the irregularities of 
equatorial electro jet has been carried out by employing the Oicamarca 
VHF radar echoes. In the present study, our interest is mainly limited 
to the case where the sporadic E occurs near the geomagnetic equator. 

The equatorial electrojet is a strong current embedded in a 
slightly ionized ionospheric plasma in the E layer near the geomagnetic 
equator. This current is a flow of electrons, motivated by E^ x B 
drift, across geomagnetic field lines through ions and neutrals, where 
E^ is the induced electric field along vertical height and B is the 
geomagnetic field. By using radar backscatter, extensive investigations 
have made it possible to elucidate physical processes for the formation 
and movement of irregularities. Balslcy (1965; 1969a, b) classified the 

*Part of the results are published in the Proceedings of International 
Conference on Recent Advances in the Physics and Chemistry of the E 
Region (Ed. by S. Matsushita and L. G. Smith) August 1974. 


irregularities into types I and II. The instabilities driven by the 
electrons drifting through ions across the magnetic field that was proposed 
by Farley (1963a, b) and Buneman (1963) are generally accepted as the 
generation mechanism of this type I irregularities. Type II irregularities 
are considered to be generated by a plasma density gradient in crossed 
electric and magnetic fields. 

The analysis of VHF radar echoes from electron concentration 
irregularities in the equatorial electrojet by Balsley (1969a) has 
indicated that there are two distinct types of irregularities in the 
region. Also, Cohen (1973), by using the spectral decomposition of the 
power spectrum of radar echoes at 50 MHz from the equatorial electrojet, 
found out that the phase velocity of type I irregularities is comparable 
to the ion-acoustic speed, and the electron drift velocity can exceed 
the ion-acoustic speed. The indication of experimental results that the 
electron drift velocity in the equatorial electrojat can be supersonic 
is contrary to some existing theoretical opinions (e.g., Sato, 1973). 

For the type II irregularities, the “xperimental results showed that 
both phase velocity of irregularities and electron drift velocity are 
smaller than the ion-acoustic speed. 

The theoretical treatment of the problem of growth of t pe I 
and type II irregularities in the electrojet has been examined extensively 
by Register (1971, 1972), Register and d'Angelo (1970), and Sudan, et al. 
(1973). However, their theoretical models were limited to the two fluid 
model (i.e., ion and electron fluid) and no neutral interaction was 
considered. Furthermore, they did not cooperate the energy transport in 
their studies and treated only the problems of the density gradients in 
their models. The problem of temperature gradients was considered by 
Cunnold (1972). However, he also used the two-fluid approach and did 
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not consider convective motions for both ions and electrons. In the present 
study, we propose an inhomogeneous three fluid model to investigate the 
unified theory of type I and II electrojet irregularities which may be 
relevant to the mechanism of equatorial sporadic E. Particularly, we are 
interested in physical mechanisms of type I and type II irregularities and 
how the electron particles are accelerated due to E x B drift in the electrojet. 
Thus, we have examined the density gradients and temperature gradients for 
three species and both as the driven mechanism for the electrojet. The 
results of computer computations based on the present model have been 
compared to the experimental observations with reasonably good agreements. 


HI-2 Mathematical Formulations 


This study is concerned with the propagation of ionospheric 
disturbances through the inhomogeneous medium which includes the species 
of electron, ion and neutral particles. Thus, a three-fluid model is 
proposed to investigate the characteristics and instabilities of the 
propagation of waves and particles. 

Let n, V, T, P and \> denote number density, velocity, temperature, 
pressure and collision frequency, respectively, and let superscripts and 
subscripts e, i and n denote electrons, ions and neutral species, 
respectively. Let E and B be the electric and magnetic fields, e the ion 
charge, m the mass, and c the jpeed of light. Then the governing equations of 
particles and fields may be written as follows: 

3n . 

J + V • (n . V J = 0 j = e, i and n (1) 

dt - J "'J 


dv 




( 2 ) 


(3) 


dv. 


-ir = - s ’’i " F '' 5> 

- "•iVie IXi - Xe’ - "’i Vin <'» 



P. V • V, = 0 
J • 


with j = e, i and n 


(5) 



where d/dt = 8/3t + (v»v). Finally, of course, the electromagnetic fields 
must satisfy Maxwell's equations. 

In the present study, the time scale of Interest is on the order 
of 0.1 second. Under this condition, momentum exchange through particle 
collisions between neutral-ion, electron-neutral, electron-ion, ion- 
neutral, and ion-electron, but not neutral -electron, plays the significant 
role in our analysis. Energy exchange through collision has bet-n ignored 
in the present model because the time scale of interest of energy exchange 
through collisions is much longer than the time scale of interest in the 
present study. 

In general, physical parameters in the upper atmosphere are in- 
homogeneous. Particularly, the inhomogeneities of number densities and 
temperatures for three species of particles along the vertical height is 
significant. As shown in Figure 1* wa choose a Cartesian coordinate 
system as a geometry for the magnetic equator, where is vertically 
upward, e^ is eastward, and e„ is northward. The inhomoqeneities of 
number densities and temperatures are assumed to be much smaller than 
their average values respectively so that Taylor series expansion with 
respect to the inhomogeneity is applicable (Mikhailovskii , 1967; Hung 
and Liu, 1971), i.e., 

"j ' ""f * ^nj • '*5> ^ 

with j = e, i and n 


( 6 ) 

(7) 
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where 
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<n 
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dz 


hi 


Jz = z, 


_L_ 


d<T.> 
— Hz 


2= Z, 


0 0 
and < > the average value and 6 the small fluctuation of physical parameters. 

The limit of small amplitude fluctuations, |6n/<n>|« 1, etc., Is 
considered for the calculation of dispersion relation. If the fluctuations 
are sinusoidal, I.e., proportional to 


exp [1 (k • X - wt)] 

equations (1) - (5) become, after substituting equations (6) - (7) and 
linearization, 


6n, 


jhk = (k + i Ki) * j = e, 1 and n 
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where £ is a unit vector in y-direction, Vihich is along geomagnetic north, 
*/ 

a! = (5/3) <T^>/m. is the square of sound speed of species j, and n. = 

J J J J 

eB/m-c is the cyclotron frequency of species j. 

J 

By using equations (8) - (12) together with Maxwell's equations, 
we can obtain the following solutions; 
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The expressions of 5 *, 5 j, hj. nj. nj. ?j. ?j. and 5 ! are very 
complicated, therefore, we will skip these results given by Hunt (1974). 

It is very interesting to point out that the neutral wind is 
mainly induced by the ion motion through ion-neutral collisions as 
shown in equation (13). If we substitute equation (14) in equation 
(13), we have 
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Equation (16) shows how neutral v/ind is accelerated by the electromagnetic 
fields through ion-neutral collisions. 

Now the electric current density, 6j, electric conductivity tensor, 
o, and dielectric tensor, K , are defined as 


and 


6 j = e <n> (6y. - dy^) 
Ka3 ^ ^ 17 


(17) 

(18) 


(19) 


By substituting equations (14), (15), (17), and (18) in equation (19), we 
have the expression of dielectric tensor KxB which characterizes the 
material constant of ionosphere, i.e.. 
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where Wpg = 4it<n> e*/nig. 

As is well known, the requirement that equation (20) be consistent 
with Maxwell's equations gives the dispersion relation 

det (k k - k“l ) + ^ ] = 0 (21) 

The dispersion relation which governs the behavior of wave-like 
disturbances In the ionosphere is quite complicated. The numerical 
computation Is generally requested to solve this problem. 
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ni-3 Preliminary Results of Numerical Computation and Conclusions 
As mentioned earlier, the purpose of the present study is to 
investigate the plasma drift instabilities in which the growing drift 
waves propagating across the field lines may contribute to the triggering 
mechanism of sporadic E in the region of magnetic equator. Since the 
dynamics of wave propagation is solely governed by the dispersion relation, 
equation (21) which is quite complicated, the numerical computation is 
a necessity to the present investigation. 

At an altitude of 110 Km, the parameters v/ith Tg= Tj^ = 250'’K, 

Vgn = 2.88 X 10** sec“S = 8.9 x 10^ sec"*, = 1.8 x 10^ sec"*, 

V = 8.9 X 10“'*sec"*, v j = 3.2 x 10"® sec"*, and Jl. = 1.5 x 10^ sec"* 

D6 m » 

are chosen. All the other values of physical parameters under non-disturbed 
conditions are chosen from U. S. Standard Atmosphere (1962, 1966). In 
order to have a better comparison with existing models of type I and type 
II irregularities, calculations are proceeded in the following three 
cases: 

(1) Case A: Electrojet Without Density Gradients 

lAr 

In this case, we assume that 0 and k-j-j f 0. 

The date of temperature gradients for three species are 
chosen from U. S. Standard Atmosphere (1962, 1966). 

Three roots with positive growth rate were obtained 
for this case from dispersion relation (eq. 21). 

Table I shows the result of this case. Westward 
drift waves with phase velocities from 303 m/sec to 
384 m/sec are excited. Growth rates of the excited 
waves (Imu/Rew) range from 0.06 to 0.17. This 
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kind of instability resembles type I Irregularities, 

(2) Case B; Electrojet Without Temperature Gradients 

★ ★ 

In this case, we assume that '»■ 0 and k^j f 0. Again 
the data of density gradients for three species are chosen 
from U. S. Standard Atmosphere (1962, 1966). Four roots 
with positive growth rate were obtained for this case from 
dispersion relation (eq. 21), Table 2 indicates this 
result. Again westward drift waves with phase velocities 
from 329 m/sec to 395 m/sec are excited. Growth rates 
of the excited waves (Imw/Rew) range from 0.04 to 0.12. 
Similar to Case A, this kind of Instability resembles 
type I irregularities. 

(3) Case C; Electrojet With Density and Temperature Gradients 
In this case, both k^j 0 and k^^j f 0. Again data for 
both temperature and density gradients for three species 
are chosen from U. S. Standard Atmosphere (1962, 1966). 

Nine roots with positive growth rate were obtained from 
dispersion relation (eq. 21) for this case. Table 3 
shows this result. Westward drift waves with phase 
velocities from 180 m/sec to 307 m/sec are excited. 

In this case, none of the v/aves, which propagate faster 
than ion acoustic speed, are excited. Growth rates of the 
excited waves are generally very high, ranging from 1 to 3. 
This instability resembles type II irregularities. 
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Comparing the results in Tables 1 to 3, some Interesting points 
are discovered. Type I irregularities of equatorial electrojet in our 
model can be generated in the inhomogeneous medium along the vertical 
height if either the temperature gradient or density gradient vanishes. 

In the meanwhile, type II irregularities can be excited if neither the 
temperature gradient nor the density gradient vanishes. Furthermore, 
the growth rate of irregularities obtained from the calculation suggests 
that the wave-wave nonlinear interaction should be taken into account for 
type II irregularities while it may be proper to proceed linearized theory 
for type I irregularities. In other words, turbulent mechanism, which is 
different from the present model, may be the characteristic of type II 
irregularities, while the laminar theory is applicable to type I 
irregularities. 

Sugiura and Poros (1969) suggested that the motion of east-west 

electron electrojet drift velocity 6v^ can be determined as 

6 

= (z). (22) 

In our expression, and In equation (15), are practically [ible 
since and are several orders smaller than Therefore, we have 

In our case, numerical computation shows that is practically 
on the order of unity which also agrees with equation (22). Since the 
vertical polarization field varies from 20 mV/m in the daytime to -20mV/m 
at night (Kato, 1973; Lee and Kennel, 1973), change of direction of electron 
drift in the equatorial electrojet from westward during daytime to east- 
ward during night hours can be fully explained based on this mechanism. 
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To check the observations made by Cohen (1973) that the 
horizontal electron velocity in the equatorial electrojet may exceed 
the ion acoustic speed, numerical computations are also made by employing 
equation (15) to calculate the westward electron velocity in the electro- 
jet for the type I irregularities. The results show that the westward 
velocity of electrons may reach 450 m/sec, which exceeds the ion acoustic 
speed, for vertical polarization electric field being 18 mV/m and west- 
ward polarization electric field being -1 mV/m. This result agrees with 
the experimental observations made by Cohen (1973). 

In closing, the following conclusions are made from the present 
study; (1) triggering mechanism for the sporadic E near magnetic 
equator is due to the inhomogeneous medium (gradients of densities and 
temperatures in three species of particles along the vertical height) 
which create cross-field plasma irregularities imposed on westward 
electrojet; (2) type I Irregularities of equatorial electrojet can be 
excited in the inhomogeneous medium along the vertical height if either 
the temperature gradient or density gradient of three species of 
particles vanishes; (3) type II irregularities of electrojet can be 
excited for the existence of both temperature and density gradients 
of three species along vertical height; (4) strong growth rate of type 
II irregularities suggests that turbulent mechanism may be the main 
structure for the type II irregularities; (5) rather weak growth rate 
in type I irregularities indicates that linearized theory may be 
applicable for type I irregularities; (6) vertical polarization field 
which shifts from positive in daytime to negative at night indicates 
horizontal electron drift in electrojet change from westward in daytime 
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to eastward at night; (7) the westward electron drift velocity In electrojet 
can be on the order of or exceeds Ion acoustic speed for type I irregularities. 

Finally, we have shown that the three-fluid model can be used to 
predict both type I and type II irregularities by Introducing proper physical 
mechanisms; namely, density gradients and temperature gradients of the three 
species. 


37 



REFERENCES 


Balsley, B. B. (1965), "Some Additional Features of Radar Returns from 
the Equatorial Electrojet," J. Geophvs. Res . 70, 3175-3182. 

Balsley, B. B. (1969a), "Some Characteristics of Non-Two-Stream 
Irregularities in the Equatorial Electrojet," J. Geophys. Res, 74. 
2333-2347. 

Balsley, B. B, (1969b), "Measurement of Electron Drift Velocities in 
the Night-Time Equatorial Electrojet," J. Atmos. Terr. Phvs . 

475-478. 

Buneman, 0. (1963), "Excitation of Field Aligned Sound Waves by Electron 
Streams," Phys. Rev. Lett. 10 , 235. 

Cohen, R. (1973), "Phase Velocities of Irregularities in the Equatorial 
Electrojet," J. Geophvs. Res. 78, 2222-2231. 

Cunnold, D. M, (1972), "An Electron Temperature Gradient Instability and 
Its Possible Application to the Ionosphere, J. Geophys. Res ., 77 , 223-233. 

Farley, D. T. (1963a), Two-Stream Plasma Instability as a Source of 
Irregularities in the Ionosphere," Phys. Rev. Lett. 10, 279. 

Farley, 0. T. (1963b), "A Plasma Instability Resulting in Field-Aligned 
Irregularities in the Ionosphere," J, Geophvs. Res. 68 , 6083. 

Hung, R. J. and V. C. Liu (1971), "Effect of Temperature Gradient on 
the Drift Instabilities in a High-$ Collisionless Plasma," Phys. Fluids 
14.2061-2063.^ 

Hung, R. J. (1974), "Dynamics of ULF Waves in the Magnetosphere," 

Research Report, University of Alabama in Huntsville. 

Kato, S. (1973), "Movement of Irregularities in the Equatorial E-Region - 
A Review," J. Atmos. Terr. Phvs . 35 , 1073-1082. 

Lee, K. , and C. F. Kennel 1 (1973), "Convective Amplification of Type I 
Irregularities in the Equatorial Electrojet," J. Geophys. Res. 78-, 4619- 
4629. 

Matsushita, S. (1951), "Intense E_ Ionization Near the Magnetic Equator," 
J. Geomaqn. Geoelect. 3, 44-46. 

Mikhailovskii , A. B. (1967), in Reviews of Plasma Physics. Vol . .3, 

(ed. by M. A. Leontovich), Consultants Bureau, New York. 

Register, Andre (1971), "Nonlinear Theory of Type I Irregularities in the 
Equatorial Electrojet," J. Geophys. Res ., 76 , 7754-7760. 


38 



Register, Andre (1972), "Nonlinear Theory of Cross Field Instability With 
Application to the Equatorial Electrojet," J. Geophys, Res ., 7£, 2975-2981, 

Register, A. and N, D'Angelo (1970), "Type II Irregularities in the 
Equatorial Electrojet, " J, Geophys. Res ., 75 , 3879-3887. 

Sato, T. (1973), "Unified Theory of Type I and II Irregularities in the 
Equatorial Electrojet." J. Geophys. Res ., 78, 2232-2243. 

Sundan, R. N., J. Akinrimisi and D. T. Farley (1973), "Generation of 
Small-Scale Irregularities in the Equatorial Electrojet," J. Geophys . 

Rer» ., 78, 240-248. 

Sugiura, M. and D. J. Poros (1969), "An Important Model Equatorial Electrojet 
With a Meridional Current System," J . Geophys. Res ., 74 , 4025-4034. 

U. S. Standard Atmosphere (1962). 

U. S. Standard Atmosphere (1966). 


39 



FIGURE CAPTION 


Figure 1 Geometry of Cartesian Coordinate System Near Geomagnetic 
Equator. 
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TABLE Ill-i CHARACTERISTICS OF ELECTROJET WITHOUT DENSITY GRADIENT 



1 



TABI/E III-2 CHAEACTERISTICS OF ELECTROJET WITHOUT TEMPERATURE GRADIENT 
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TABLE III-3 CHARACTERISTICS OF ELECTROJET WITH DENSITY AND TEMPERATURE GRADIENTS 





APPENDIX 

PARTICLE VELOCITIES AND MOBILITY TENSORS IN 
INHOMOGENEOUS THREE-FLUID MODEL 


Dynamics of particle acceleration and wave excitation are governed 
by the fundamental equations of three-fluid model shown in equations (1)- 

(5) . If we assume that the inhomogeneities of temperatures and densities 
for three species of particles are mainly along the vertical height, we 
have the expressions equations (6) and (7). If we consider only the 
small amplitude fluctuations which are sinusoidal, fundamental equations 
can be rewritten as equations (8) - (12) after substituting* equations 

(6) and (7) and linearization. To solve the equations (8) - (12) together 
with Maxwell's equations, let us assume that the waves are mainly propa- 
gating in north-east plane. Then, the velocity of neutral particles can 
be obtained in terms of the velocity of ion particles 
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Here, subscripts and || denote the components perpendicular and parallel 
to the geomagnetic field, respectively. In the present model, subscripts 
_L and [I particularly imply the components along east and north directions, 
respectively. 

Velocity of ion particles can be solved in terms of mobility tensor, 

i.e.. 
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where 
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Finally, velocity of electron particles also can be solved in terms 
of mobility tensor, i .e. , 
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CHAPTER IV 


TIME-DEPENDENT NEUTRAL WIND STRUCTURE 
IN THE AURORAL E REGION* 


IV-1 Introduction 

Recently, Brekke.etal. (1973) presented neutral winds In the 
auroral E Region, deduced from measurements of the ion drift velocity 
at different altitudes by the incoherent scatter radar facility at 
Chatanika, Alaska. Because of experimental limitations, only two 
data points per profile, corresponding to 110 km and 167 km, lie 
approximately within the E Region. TKese values are weighted averages 
over several tens of kilometers in altitude so that no structure between 
these altitudes is available. In this study a theoretical model is 
constructed in an attempt to determine the neutral wind structure in 
this region. 


*This work was done by S. T. Wu and R. H. Comfort. Portions of the 
results are published in the Proceedings of International Conference 
on Recent Advances in the Physics and Chemistry of the E-Region (Ed. 
by S. Matsushita and L. G. Smith), August 1974. 
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IV- 2 Analysis 


The model developed Is based on electric fields as the primary 
driving force for the Ions, with momentum transferred from the Ions to 
the neutral gas through collisions driving the neutral wind. Both Redder 
and Banks (1972) and Heaps (1972) have demonstrated that viscosity In the 
neutral gas Is important, even at E Region altitudes, In redistributing 
momentum vertically; viscosity is therefore Included In the neutral 
equation of motion. In addition, Heaps (1972) has shown that the 
Coriolis force Is Important in changing the direction of the neutral 
wind over the time scale of Interest (several hours); this force is 
also Included in the neutral equation of motion. In this initial 
state of model development pressure gradient terms are neglected, as 
are the continuity and energy equations. Model values are used for 
the densities and temperatures governed by these equations. Gravi- 
tational forces. are assumed to be balanced by vertical pressure gradients 
In the static diffusion models from which model values are obtained; 
therefore, as a first approximation, gravity and the vertical pressure 
gradient are not included explicitly in the vertical equations of motion. 

With these assumptions, the governing equations for the present 

m 

problem are the ion equation of motion 
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and the neutral equation of motion 
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where t is the stress tensor, given by 
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In equations (1) - (3) and are the ion and neutral wind velocities 
Vin and are average ion-neutral and neutral -ion collision frequencies 

E and B are the electric and geomagnetic field strengths, and are 
the mean ion mass and neutral mass density, n is the coefficient of 
absolute viscosity, e and c are the electronic charge (magnitude) and 
speed of light and is the earth's angular velocity vector. Horizontal 
gradients in these quantities are assumed to be negligibly small. 

Average collision frequencies and are evaluated in the following 
way: 

^in " S ^r '^rs ^rs 

- ■ Pi - 
''ni " ^'’in 


where r and s range over all ion and neutral species respectively, p 

I 3 

is the reduced mass, and is the von-neutral collision frequency for 
momentum transfer as given by Banks (1966). 

For convenience in performing the calculations, the ion equation 
is treated in the "plasma" coordinate system (Figure la) and the neutral 
equation in a local geographic coordinate system (Figure lb). The 
transformations between these coordinate systems for an arbitrary vector 
q are: 

(i) From plasma coordinates (1, 2, 3) to local geographic 
coordinates (x, y, z) 
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(4a) 

(4b) 

(4c) 


« qg s1n6 - cos6 (q^ sin I + q^ cos I) 

qy = qg cos6 + sin6 (q-j sin I + q^ cos I) 

(11) From local geographic coordinates (x, y, z) to plasma 
coordinates (1, 2, 3) 
q-| = q^ cos I - sin I (q^^ cos6 - qy s1n6) (5a) 

^2 ° % s1n6 + qy cos6 (5b) 

qg = -q^, sin I - cos I (q^^ cos6 - qy s1n6) (5c) 

The angles I (magnetic dip angle) and 6 (magnetic declination angle) 
are defined by Figure Ic. All results are presented In the local 
geographic coordinate system. 

Because of the high conductivity along magnetic field lines, 
the electric field t Is assumed to have no component parallel to the 
magnetic fields and E Is assumed to be mapped down the field lines 
without loss or rotation to all altitudes of interest. Brekke, et al., 
(1973) provides data for the Ej^ and Ey components of electric field. 

The assumption of no parallel electric field (E^ = 0) allows equation 
(5c) to be solved for E^, giving 

= - cot I (Ej^ cos6 - Ey s1n6) . (6) 

Since E Is used only in the plasma coordinate system (1, 2, 3), equation 
(6) Is used together with equation (5) to express E-j and Eg In terms 
of the E and E data. 

A J 

In early numerical tests, using an explicit finite difference 
technique. It was determined that the Ion equations of motion come 
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into a steady state with constant driving forces within a few collision 
times models used, these collision times vary from 

-10"'' sec at 90 km to -1 sec at 250 km. Hence, at all altitudes of 
Interest the Ions achieve a steady state with driving forces within a 
few seconds, providing the driving forces vary only slightly during 
this time. In this problem, the driving force is an electric field 
which varies significantly only over a period of several minutes (values 
presented by Brekke, et al. (1973) represent at least 10 minute averages). 
The time-dependent Ion equations can therefore be replaced by the time- 


independent equations 


(■ 


at 


« O) with little loss of accuracy. Ion 


component equations in plasma coordinates then have the customary form 





1 + 
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(7a) 


(7b) 


''13 ' %3 

where = (eB)/(fnjc) is the mean Ion gyrofrequency. The time dependence 

of v^ is now implicit in the time dependence of E and v^^. 

The coupled set of time-'independent ion momentum equations (7a)- 

(7c) and neutral momentum equations (2) is solved numerically, using an 

explicit finite difference technique, by prescribing a value of the 

electric field at each time step. As initial conditions, all neutral 

velocity components are taken to be zero. Boundary conditions at 90 km 


are the vanishing of all neutral velocity components, and at 250 km, 

3v 

>r-^=0, These boundary conditions, while not strictly correct, have 

oZ 

little effect in the region of interest, 100 km to 200 km. Two different 
sets of electric field values are considered, a constant electric field 
and electric field values inferred by Brekke, et al., from data taken 
with the incoherent scatter radar at Chatanika, Alaska, Model neutral 
parameters used are obtained from a Jacchia-71 model atmosphere at a 
constant exospheric temperature of 1000 degrees. Ion densit-!es, both 
day and night values, are taken from Brekke, et al. (1973). Crude 
ion composition models for day and night have been formed based on the 
rocket data of Kopp, et al, (1973). The only time variation associated 
with these models in the calculations is the use of the night model 
during hours of darkness for the calculation using Chatanika data, 
while using the daytime model for all other times. With these models, 
at 122 km. 

IV-3 Results 

In order to understand more clearly how the ion-neutral system 
responds to electric fields, a simple case is considered initially. A 
co.,stant electric field of 20 mV/m is applied for five hours in the 
+ X (south) direction, followed by three hours of zero field. All 
geometric parameters correspond to the location of Chatanika, Alaska 
(next example), and daytime model values are used. Figure 2 shows 
the ion and neutral velocity components as functions of time at 
altitudes of 115 km, 125 km, and 150 km. The difference in magnitude 
of the ion components at these altitudes clearly displays the effects 
of increased colli si onal coupling with the neutral gas at lower 
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altitudes. In addition, tho Inertia of the neutral gas Is evident at 
T = 5 hours, when the electric field driving force Is terminated; the 
"coasting" neutral wind then becomes the driving force for the ions. 

To see the altitude structure of the velocity fields, altitude 
profiles of all velocity components are shown for T = 5 hours In Figure 
3a. The same calculation has been performed with the corlolls term 
set to zero in order to determine the effect of that term on the calcu- 
lation; resulting profiles are shown In Figure 3b. Comparison with 
Figure 3a shows clearly that the corlolls force Is primarily responsible 
for the magnitude of the vertical components. This force Is also seen 
to cause a redistribution of momentum among the horizontal neutral velocity 
components, with consequent effects on the Ion velocity. A similar 
calculation setting the viscous terms to zero gives results (not shown) 
very similar to Figure 3a. The primary effects of viscous terms are 
found to be a small damping of the vertical components and a minor 
smoothing of the altitude profiles. 

For more realistic time variations, electric fields determined 
by Brekke, et al. (1973), from Incoherent scatter radar observations 
at Chatanika, Alaska (L 5.7, A = 65”, I = 76.5® , 6 = 29®) for July 10-11, 

1972, are used. These electric field values are shown In Figure 4a (note 
that here +x1s south, while Brekke, et al. (1973), used positive 
values for north). Detenninatlon of electric field values was made 
from high altitude Ion drift measurements, while lower altitude ion 
drift measurements were used to deduce neutral winds. The technique 
for the latter determination, described in detail by Brekke, et al. 

(1973), results in a weighted height- Integrated value for neutral winds 


with the primary weighting in the 110 km to 120 km region, but varying 
with the ion density profile. These neutral wind values froffi Brekke, 
et al, (1973) are shown in Figure 4b. 

Time variation of ion and neutral wind components, calculated 
as described above for the electric field of Figure 4a, ir, shown in 
Figure 5 for altitudes of 115 km» 125 km, and 150 km. Ion velocity 
temporal structure parallels that of the electric field while the 
neutral components vary smoothly, almost sinusoidally in time. Altitude 
effects are particularly noticeable in the y components. Representative 
examples of the altitude structure are shown In Figure 6. All ion and 
neutral velocity components are shown as functions of altitude for 
0800 hours UT in Figure 6a and for 1030 hours UT in Figure 6b. These 
profiles, as those in Figure 3, are characterized by large gradients in 
the 110 km to 130 km region, and in some cases, a reversal of direction 

there. Above 140 km the velocity components change rather slowly with 

* 

altitude, reflecting the reduced collisional coupling in this region. 

A comparison of the neutral velocity components in Figure 4b 

with those in Figure 5, while not strictly appropriate quantitatively, 

• 

is nevertheless instructive. One notable difference is the relative 
lack of temporal structure In the calculated components. This 
qualitative discrepancy should have a physical explanation. If the 
calculated values for 115 km are used in the comparison, the peak 
magnitudes agree and the general directions of the x components agree; 
however, agreement between the y components is rather poor. If a 
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higher altitude is picked for the calculated y component, agreement in 
general direction is improved - at the expense, however, of agreement 
in general magnitude. 

With this comparison as a background, several shortcomings of 
these calculations can be examined qualitatively. In the present 
treatment only two representative ion density altitude profiles are 
used (for day and night conditions), whereas the data analysis of 
Brekke, et al. (1973) employed profiles measured simultaneously with 
the ion drift measurements. The high variability of these profiles 
is Illustrated in data presented by Banks, et al. (1974). Although 
inclusion of this variability (through v^^) in the neutral equation 
of motion could hvjdify the calculated neutral winds somewhat, the 
result would probably be similar to that due to electric field 
variations, which enter the neutral equations through ion velocity 
variations in the collision term. Figure S demonstrates that such 
short period variations are effectively filtered out of the neutral 
components, due to the large inertia of the neutral gas relative to 
the ion gas. Furthermore, the tendency of electric fields and E region 
ionization densities to be anti correlated (Banks, et al., 1974) would 
have a compensating effect in the collision terms of the neutral* 
equations, reducing the effective magnitude of variations in these terms. 

Another omission of the present calculation is horizontal pressure 
gradients. It is reasonable to expect large short term variations in 
such pressure gradients associated with corresponding large variations 
in Joule heating in the auroral oval. However, no data is available 
to provide information on these pressure gradients for the time and 


locjition treated. There is also little theoretical work to suggest how 
far from the source region significant pressure gradients might extend. 

In addition to the large variable pressure gradients associated with 
Joule heating, smaller pressure gradients due to diurnal thermal tides 
would bias the results, particularly when electric fields are small. 

For the present calculation, this effect is not expected to modify 
the results to a large extent. 

Based on the comparison of Figure 4b and 5 above, it appears 
that the effective altitude of the measurement may be varying. Since 
it is the ion drift which is measured, the altitude structure of the 
ion velocity components must enter into the weighting which determines 
the altitude region contributing most significantly to the measured 
results. Figures 3 and 6 show that below 110 km velocities are small, 
increasing rapidly with altitude to about 130 km and less rapidly above 
that. In addition reversals of velocity component directions in the 
no km to 130 \m region would effectively cause cancellations in the lower 
altitude contributions to the integrated results, further raising the 
effective contributing altitude region. If these shifts occur within the 
region of large vertical gradients, it appears that significant differ- 
ences in the weighted integrated results could occur. Examination of 
Figure 5 Indicates that ion velocity altitude structure can change very 
rapidly in response tc changing electric fields. This, variability 
coupled with the variability of ion density altitude profiles, both 
of which affect the altitude weighting of the integrated velocity 
measurement, could be factors in explaining the temporal structure of 
the neutral velocities deduced from measured ion velocities (Figure 
4b). These qualitative arguments require quantitative support; 
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calculations of this nature are in progress. 

Finally, the vertical velocity components of lOu m/s or more, 
evident in all of the calculated results (Figures 3, 5 and 6), require 
comment. Comparison of Figures 3a and 3b has established that the 
Coriolis term is a primary source of these velocities. However, 
measured ion vertical velocities are typically less than DO m/s (Banks, 
et al. 1974; Rieger, 1974). The assumption of total cancellation 
between the vertical pressure gradient and gravitational acceleration 
is clearly not a good one for the dynamic case. Order of magnitude 
estimates indicate that in the ion equation, these terms are smaller 
than those included by factors of 10^ to 10®. However, in the neutral 
equation, the gravity term is 10® times larger than the terms included, 
so that even a small perturbation of the hydrostatic pressure can have 
a very significant affect on vertical velocities. Pressure perturbations 
due to vertical driving forces must therefore be included in a realistic 
treatment of vertical winds; and this requires inclusion of the con- 
tinuity and energy equations in the set of equations to be solved. The 
qualitatively significant result remains, however; large horizontal 
winds mean significant vertical driving forces are present. 

IV-4 Conclusions 

Over time scales of several hours, the coriolis force plays a 
significant role in redirecting the neutral wind, including driving 
a vertical component; for these time scales this force should be 
Included in the neutral equations of motion. Ion velocities respond 
rapidly to changing electric fields at all altitudes so that vertical 
structure of these velocities can change quickly. Resulti.':g changes 
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in the neutral velocity altitude structure* however, require one to 
two hours, due to the greater neutral gas density. In these calculations, 
ion and neutral velocities are small (below 110 increase in magnitude 
rapidly with altitude to about 130 km, and change more slowly with 
altitude above. This variability of Ion velocity altitude structure 
may complicate the interpretation of weighted height-integrated ion 
velocity measurements. Realistic treatment of vertical velocities 
requires explicit consideration of pressure perturbations resulting 
from vertical motion; in general this will require inclusion of the 
continuity and energy equations in the system of equations. 

Finally, we may conclude from the present calculation that 
the effects of coriolis force are clearly indicated. A strong shear 
layer appeared in the neutral wind profiles. However, the Richardson 
number may still very well be in the laminar flow region. 
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FIGURE CAPTIONS 


Figure IV-1 Definitions of coordinate systems used in the calculation 
and the angles relating them. 

Figure IV-2 Ion and neutral velocity components as functions of time 
for altitudes 115 km, 125 km and 150 km for the constant 
electric field: E = 20 mV/m, T <. 5 hours; E = 0, 

A A 

T > 5 hours; E = 0, all T. 

(j 

Figure IV-3 Altitude profiles of ion and neutral velocity components 
at T = 5 hours for the electric field defined for Fig. 2: 
(a) all terms in equation (2); (b) omitting coriolis 
term from equation (2), 


Figure IV-4 Data determined from observations by the incoherent scatter 
radar at Chatanika, Alaska, on 10-11 July 1972 {from Brekke, 
et a1. 1973). (From 9730 to 0900 hours UT, no data was 
reported; electric fields for that period are sketched in 
graphically for continuity in the calculation.) (a) electri 
fields; (b) neutral wind velocities. 


Figure IV-5 Calculation ion and neutral velocity components as functions 
of time for altitudes 115 km, 125 km, and 150 km due to 
electric fields of Fig. 4a, 


Figure Iv-6 


Representative calculated altitude profiles of ion and 
neutral velocity components due to electric fields of Fig. 
4a. (a) for UT = 0800 hours; (b) for UT = 1030 hours. 
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Figure IV-2{a) (x-component) 
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Figure IV-2(b) (y-component) 
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Figure IV-2(c) (z-component) 
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CHAPTER V 


CONCLUDING REMARKS 

In this study we have presented a multifluids plasma model for 
the ionosphere in the region of 90-2,000 Km. The wide range of density 
and temperature distributions will lead to a wide-range distribution of 
the collision frequencies among the species (i.e., electrons, ions and 
neutrals) at different altitudes. Thus, the physics of the dynamical 
behavior in this region can be classified based on these parameters 
(i.e., density, temperature and collision frequency). Subsequently, 
the proper equations to govern the physical phenomena can be obtained. 

In the present study, we have presented only three examples 
demonstrating the theory developed which will enable us to study the 
wide range dynamical phenomena in the upper atmosphere (90- 2,000 Km), 
including the lower part of the magnetosphere. Again, there are many 
Interesting problems observed in the ionosphere and magnetosphere 
regions. We shaM mention a few which form a basis for future study. 
These are; 

i) Utilize the present model to study the joule heating 
effects in relation to the upper atmospheric dynamics, 

ii) The detailed examination of the instability mechanism 
as an altitude dependence should be very interesting; 
it could be possible to shed some light concerning 
the triggering mechanism of the ionospheric irregularities, 

iii) All the examples presented are limited to the local 

coordinate system. It is important to extend the present 
model to the global scale. 
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